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We study the critical charge dynamics of the superconducting to the normal-state transition for 
La2-2!Sr2,Cu04 (LSCO) thin films with a wide range of the Sr concentration, by measuring the 
frequency-dependent excess parts of the complex microwave conductivity, which is induced by the 
superconducting fluctuations. We present a dynamic scaling analysis of the complex fluctuation 
conductivity, which includes the information on the universality class and the dimensionality of the 
critical charge dynamics as a function of the Sr concentration, the film thickness and the magnetic 
field. In our previous study (H. Kitano et al, Phys. Rev. B 73, 092504 (2006).), the 2T>-XY critical 
dynamics for underdoped LSCO and the 3D-Xy critical dynamics for optimally doped LSCO were 
reported. In this study, we observed a novel two-dimensional unknown critical charge dynamics for 
overdoped thin films from a;=0.17 to 0.20, which is clearly distinguished from the 2D-XY critical 
dynamics. Through the systematic measurements by changing the film thickness or by applying 
small magnetic field, it was confirmed that this unusual behavior, which is referred as 2D-"U" 
below, was not induced by the finite size effect but was intrinsic to the overdoped LSCO, Thus, it 
was found that the critical behavior in the phase diagram of LSCO is classified into the following 
three types; (i) 2D-XY for underdoped region, (ii) 3D-XY for optimally doped region, and (iii) 
2D-"U" for overdoped region. In other words, the dimensionality in the critical charge dynamics 
is changed twice with hole doping. We discuss possible origins of such anomalous dimensional 
crossovers with hole doping, including an interpretation based on the possible existence of a hidden 
quantum critical point near the optimally doped region. 



I. INTRODUCTION 

One of the most striking features in the cuprate super- 
conductors is the strong doping dependence of the critical 
temperature, Tc- It is well known that a plot of Tc versus 
the carrier doping forms a bell-shaped phase diagram. In 
spite of a large number of experimental and theoretical 
studies on this issue, the physical origin of the phase di- 
agram including the pseudogap phenomena above Tc in 
the underdoped region is still debated. 

When we focus on a role of quantum criticality (QC) in 
the phase diagram of high- Tc cuprates, most of theoreti- 
cally proposed phase diagrams can be classified into the 
following two groups. In the first group, two quantum 
critical points (QCPs) arc recognized at both endpoints 
of the superconducting transition line. However, the bell- 
shaped superconducting dome or anomalous properties in 
the pseudogap regime is explained by a different concept 



from the QC, for instance, the gauge field fluctuations 
in the t — J model, ^'^ or the classical phase fluctuations 
of superconducting orders, which are assumed to survive 
even in the pseudogap regime 

On the other hand, in the second group, in addition 
to the two well established QCPs in the first group, an- 
other QCP is assumed to exist (or to be hidden) inside 
the superconducting dome. In other words, the exis- 
tence of a hidden order is assumed, which may compete 
or cooperate with the superconducting order. In this 
view point, superconductivity and the anomalous behav- 
ior in the pseudogap regime are related to the hidden 
order and its large quantum fluctuation around the hid- 
den QCP. Various candidates for this hidden order have 
been proposed; for instance, a magnetic Neel order^^ 
an incomensurate charge-density-wave;^ a time-reversal- 
violating state,— a d-density wave,— and a charge/spin 
stripe order Possible candidates for such an unknown 
order have also been investigated by a more general con- 
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sideration using a group-theoretic classification for com- 
peting orders. Several experiments have suggested a 
possibility that there is the hidden QCP in the vicin- 
ity of the optimally doped concentration. ^^si^iHii^ii^iii 
However, the interpretation of such experimental results 
is still controversial, and the relationship with a mecha- 
nism of superconductivity also remains unresolved. 

These two groups contrast strikingly with each other 
in term of the possibility of a hidden quantum phase 
transition (QPT). If a typical frequency scale of such 
a QPT was sufficiently larger than a thermal energy, 
it is difhcult to disentangle quantum and thermal ef- 
fects This means that the classical critical dynam- 
ics can be affected by the quantum-fluctuations effects 
near the hidden QCP. Therefore, one can know which 
group is more appropriate for describing the phase di- 
agram of high- Tc cuprates, by investigating the super- 
conducting to the normal-state transition with a wide 
range of carrier doping. Fortunately, the critical fluctu- 
ation effects, which reflect the transition nature, can be 
explored in the high- Tc cuprates, because the short coher- 
ence length, the small superfluid density and the quasi- 
two-dimensionality largely enhance the classical fluctua- 
tions of the superconducting order 

In our previous paper^ we reported that the critical 
charge dynamics for underdoped La2_a;Sra;Cu04 (LSCO) 
from x=0.07 to 0.14 was successfully expressed by the 
2D-XY universality class, that is, Berezinskii-Kosterlitz- 
Thouless (BKT) picture. On the other hand, the criti- 
cal behaviors for almost optimally doped LSCO (a;=0.16) 
were found to be expressed by the 3D-XY universality 
class. This sudden change of the universality class was 
highly unexpected, because the classical critical behavior 
must be universal, irrespective to the microscopic details, 
such as a carrier concentration.^^ 

These anomalous results supported the idea that the 
systematic measurements of the critical charge dynamics 
as functions of carrier concentration are quite effective to 
study the phase diagram of high- Tc cuprates. However, 
the previous work covered a part of the whole phase di- 
agram. Thus, it is crucially important to complete the 
investigation including the overdoped region. 

Based on these backgrounds, we have investigated the 
critical charge dynamics of the superconducting to the 
normal-state transition for high-quality LSCO thin films 
with a wide range of the Sr concentration. We measured 
the frequency dependence of the complex microwave con- 
ductivity, cr[uj) — (Ji{lu) — 1(72 (w), which is enhanced by 
the superconducting fluctuation, under zero and small 
finite magnetic fields. In this paper, we present new re- 
sults of the critical charge dynamics for optimally doped 
LSCO (a;=0.15 and 0.16) and overdoped LSCO (from 
a;=0.17 to 0.20). The most important finding was that 
there is another dimensional crossover from 3D to 2D be- 
tween a;=0.16 and a:=0.18. To clarify the details of this 
anomalous behavior, the finite size effects were also in- 
vestigated, by changing the film thickness or by applying 
small finite magnetic fields. Together with our previous 



study,^° we discuss a whole picture of the critical charge 
dynamics in the phase diagram of LSCO in terms of var- 
ious proposed models. We consider that the newly found 
second sudden crossover needs to be explained by assum- 
ing the additional hidden QCP near the optimal doping. 

This paper is organized as follows. In Sec. II, we sum- 
marize the practical merits of a dynamic scaling analy- 
sis using the frequency-dependent complex conductivity, 
comparing with the scaling analyses using other thermo- 
dynamic and transport properties. In Sec. HI, we briefly 
describe the preparation of high-quality LSCO thin fllms 
and the broadband technique to obtain the frequency- 
dependent microwave conductivity. Section IV contains 
all the experimental results on the dependence of the crit- 
ical charge dynamics on the hole-doping, the film thick- 
ness and the magnetic field. In Sec. V, various proposed 
models are systematically discussed to explain the sets of 
our data. Finally, we conclude this work in Sec. VI. 



II. DYNAMIC SCALING ANALYSIS OF 
COMPLEX FLUCTUATION CONDUCTIVITY 

The critical behaviors are characterized by the diver- 
gence of typical length and time scales in the vicinity of a 
critical point. In the case of classical phase transition^^^ 
a correlation length, ^, is often used as the typical length 
scale. ^ diverges as ^ = £,o\T/Tc — where £_o is the 

correlation length at T = and v is a static critical expo- 
nent. On the other hand, the divergence of a correlation 
time, T, is given by r = Here, z is a dynamic critical 
exponent. The divergence of ^ and t in quantum phase 
transition can also be defined in a similar manner, except 
that a singular point is not a Tc but a QCP at T = 0. 

Fisher, Fisher, and Huse (FFH)^^ provided a general 
formulation of the dynamic scaling hypothesis for the 
frequency-dependent complex fluctuation conductivity, 
CTfl(w), near the superconducting transition, as follows, 

afl(c.) « r+'"'5K^). (1) 

Here, S{x) is a complex universal scaling function and d 
is a spatial dimension. As was emphasized in our pre- 
vious paperf2£ the most essential part of our analyses 
is that we can check the validity of this hypothesis, by 
measuring the frequency dependence of (Tfl(a;). Follow- 
ing a pioneering work by Booth et al.,^^ we use both the 
magnitude, \(t\, and the phase, 4>a{= tan^^[fT2/o'i])i of 
afi{uj) as scaled quantities in the scaling analysis of (Tfl(tj). 
First of all, the phase part of afi{uj) at different temper- 
atures is scaled by using a normalizing factor, uiq, along 
the cj-direction in the plots of (j^^r versus lu. Next, the 
magnitude part of an as a function of the normalized fre- 
quency, Lij/ujQ, is also scaled by using another normalizing 
factor, ao, along the |(T|-direction. Note that ujq and ctq 
are independently obtained in our procedure, since we 
measure data sets of crfi(w) as complex quantities with 
two independent components. 
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If the data sets of and those of |(T|/cro collapse on to 
a single universal function of lo/liJo, respectively, we find 
that the dynamic scaling hypothesis is satisfied. Thus, 
Eq. Il]) tells us that the two scaling parameters, ujq and 
(To, can be connected with the diverging length scale, ^, 
as follows, 

(X « |T/r, - (2) 
(TO cx oc \TlTc - (3) 

By using these relationships, we can directly determine 
the two critical exponents, v and z, and the dimension- 
ality, c?, through the temperature dependence of ojq and 
(To. On the other hand, if the data sets of 4>a or those 
of |(t|/(To do not collapse on to a single curve, we find 
that the dynamic scaling hypothesis is broken down, so 
that we cannot extract the information of the universal- 
ity class from Eqs. (2) and (3). Indeed, we have encoun- 
tered such a breakdown of the data collapse in 4>„ for the 
conventional NbN superconducting thick film, showing a 
dimensional crossover due to the divergence of ^ beyond 
the film thicknessi^ We also observed that a broad distri- 
bution of Tc caused a breakdown of the dynamic scaling, 
which will be described in Sec. V. 

The dynamic scaling procedures described above yield 
more reliable information than other scaling procedures, 
because of the following reasons. First of all, the scal- 
ing procedure of 4)a{^) was performed in a linear scale. 
This required a more strict criterion in data collapse 
than ordinary scaling procedures performed in the log- 
log plots, which easily lead to spurious scaled behaviors. 
In addition, the dynamic scaling analyses can be per- 
formed by using only the experimental data. In other 
words, we can check whether the data obey the scal- 
ing hypothesis or not, without any assumption. This is 
in contrast with other scaling procedures. For example, 
in the case of the scaling analysis of the I-V character- 
isticSf^ the dimensionality, d, has to be assumed first. 
The value of z is obtained from the limiting behavior of 
V{(x at Tc. Only after those procedures, the 

data can be scaled in the plots of {I/T)\T - Tc\-'^'^d-i} 
versus {V/ I)\T -Tc\-'''^'^+''-''\ The value of i/ is obtained 
from a successfully scaled behavior of the plots. However, 
the ambiguity of Tc easily gives rise to large error bars on 
the obtained results of z and i/, which degrades the re- 
liability of the scaling analysis. In addition, one cannot 
detect the breakdown of the dynamic scaling hypothe- 
sis due to the dimensional crossover in these procedures. 
Other scaling analyses require the assumpition of a par- 
ticular model in advance to perform the scaling analyses 
as well, such as, the scaling analysis of the dc magnetiza- 
tion;^ the specific heat^ the thermal expansion^ and 
the ac conductivity at a fixed frequencyjS£ 

Unfortunately, there was no consensus among the pre- 
vious results on the superconducting fluctuations of high- 
Tc cuprates, mostly based on the measurements described 
above. We believe that such controversial results were 
attributed to the lack of the experimental check on the 
applicability of the static or dynamic scaling hypothesis. 



In this sense, only the dynamic scaling analysis of (ta {uj) 
is exceptional, since no assumption is needed to check 
whether the data sets oi (j)^ or \a\ collapse on to a single 
curve or not. This is the reason why we adopted this 
method to investigate the critical charge dynamics near 
the superconducting transition of high- Tc cuprates as a 
function of carrier doping. 

More precisely speaking, the determination of v and z 
is affected by the ambiguity of even if the dynamic 
scaling form predicted by FFH is confirmed to be satis- 
fied. However, the determination of the dimensionality, 
d, is free from the ambiguity of since a product of coq 
and (To is proportional to as derived from Eqs. (2) 

and (3). This suggests that LUoao is independent of T for 
d=2, while it increases with increasing T for d=3, since 
1/ is always positive. In addition, Eq. ([1]) shows us that 
the phase angle and the magnitude of (Tfl(w) behave as. 



respectively, as T approaches T^. This limiting behav- 
ior is also very useful to discuss the dimensionality. If 
d=2, then ~ tt/2 and \<t\ oc 1/w for T — s- Tc, inde- 
pendent of z. On the other hand, for d—3, it is expected 
that (t>a ~ 7r/4 and \a\ oc l/w^ for T Tc, since a fully 
relaxational dynamics with z«2 (that is, the so-called 
model-^ dynamics in the Hohenberg and Halperin clas- 
sification^) is expected for a charged superfluidfiS Thus, 
we can discuss the dimensionality of the superconduct- 
ing phase transition through the limiting behavior of (j)^ 
and loqcto in the vicinity of Tc, even if the precise value of 
Tc is unknown. This is particularly important when we 
discuss the dimensionality as a function of carrier doping. 

III. EXPERIMENTAL 
A. High quality LSCO films 

Epitaxial LSCO thin films covering a wide range from 
a::=0.07 to 0.20 were grown on LaSrA104 (001) substrates 
by a pulsed laser deposition technique. Details of the 
growth condition were described elsewhere i^i^i^ Pure 
ozone in an atmosphere of 10 mPa was used for oxidation 
during both the deposition and postannealing processes 
in the growth chamber, in order to obtain superconduct- 
ing LSCO films with very low residual resistivity. The 
excellent c-axis orientation of each film was confirmed by 
the X-ray diffraction. The rocking curve of the 002 re- 
flection showed the full width at half maximum (FWHM) 
of about 0.2°, which was almost the same as FWHM of 
LaSrA104 (LSAO) substrate, suggesting that the mosaic- 
ness of the film is not so much different from that of the 
substrate. '^^ 

All films were carefully annealed in air (or O2 gas for 
x=0.20) at 700 °C for 2-4 hours, following the heat treat- 
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TABLE L Various parameters of LSCO films. See the text 
for details. 



X 


t (nm) 


p (mficrn) 


T« (K) 




^scale ^ 


0.07 


460 


0.78 


18.0 


3.6 


18.2 


0.10 


220 


0.35 


28.2 


1.7 


28.5 


0.12 


230 


0.28 


33.7 


1.5 


33.7 


0.14 


270 


0.19 


38.9 


0.5 


38.9 


0.15 


150 


0.22 


31.6 


1.8 


32.1 


0.16 


140 


0.12 


35.5 


1.6 


35.9 


0.17 


140 


0.15 


34.4 


1.7 


34.7 


0.17 


115 


0.14 


34.4 


1.2 


34.7 


0.18 


240 


0.11 


34.5 


1.3 


34.8 


0.18 


120 


0.15 


27.9 


2.7 




0.18 


60 


0.15 


27.8 


1.5 


28.2 


0.19 


135 


0.09 


26.7 


1.0 


26.9 


0.20 


90 


0.05 


32.4 


1.7 


32.8 



ment procedure developed for high-quality LSCO single 
crystals. '^^ The annealed samples were rapidly quenched 
down to low temperatures using liquid nitrogen. 

The in-plane dc resistivity, pdc, of each film was mea- 
sured using a standard four-probe method. Table I lists 
a thickness, the resistivity at T=50 K, Tc defined by zero 
resistance, ATc defined by the 10%-90% transition width, 
and Tc used in the dynamic scaling analysis for all the 
samples. As shown in Fig. [1] and Table I, we confirmed 
that the absolute value of pdc just above Tc systemat- 
ically decreased with increasing hole doping within the 
experimental error, and also that it agreed with the re- 
ported best values for LSCO thin films^^''^^ and single 
crystals^-i with the same value of x, within the factor 
of 2. Together with the narrow rocking curve observed 
by the x-ray diffraction, these properties show that the 
thin films under this study are of so high quality that the 
critical phenomenon near Tc deserved to be investigated 
experimentally. Note that the transitions appear to be 
broad, especially for underdoped films. However, this is 
a consequence of the fluctuation, and is not due to the 
inhomogeneity of Tc- The only exception is the 120 nm 
thick film with x—0.18, which will be referred in Sec. V 
as an example of the breakdown of the dynamics scaling 
due to the inhomogeneity of Tc- 

As was already described in our previous paper for 
the purpose of investigating the critical charge dynamics 
in a whole region of the phase diagram, LSCO is an ideal 
system with a simple layered structure, where the hole 
concentration can be widely controlled. In addition, the 
use of LS AO substrate for LSCO films provides an advan- 
tage over bulk crystals of LSCO, because the tetragonal 
symmetry of LSAO substrate prevents the corrugation in 
the Cu02 plane, which is attributed to a staggered rota- 
tion of CuOg octahedra in low-temperature orthorhom- 
bic (LTO) phase of bulk LSCO.^^ Thus, we can expect 
that ideal fiat square Cu02 planes, which are free from 
disorders due to corrugations and twin boundaries, are 
realized in LSCO films on LSAO substrate. 



B. Microwave broadband technique 

When the film thickness, t, is sufficiently smaller than 
the skin depth, S, the frequency-dependent microwave 
complex conductivity, cr(Li;), can be obtained from the 
complex reflection coefficient, Sii{uj), as follows, 
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(6) 



where t and Zo{= 377 fl) are the film thickness and the 
impedance of free space, respectively. 

We measured the frequency dependence of 511(0;) in 
the frequency range between 45 MHz and 10 GHz, by 
using a vector network analyzer (HP 8510C) in a step- 
sweep mode. The thin film sample with the electrodes 
in a Corbino disk shape was terminated to the end of 
a coaxial cable through a modified jack-to-jack coaxial 
adaptor. A stable electrical contact between the sample 
and the modified coaxial adaptor was realized by two 
springs, following the work by Booth et alM' One spring 
was embedded into a center conductor pin of the adaptor, 
and the other spring was located at the backside of the 
LSAO substrate Other details of the experimental setup 
will be described elsewhere.— 

In practice, the measured reflection coefficient, 
S^^^^(^uj), is affected by the attenuation and the phase 
shifts in the intervening coaxial cable, giving rise to large 
systematic errors at lower temperatures and higher fre- 
quencies. Thus, S'J°°'^''(a;) is given by the following equa- 
tion,ii 



C*mc 
"^11 



Eji{uj)Sn{uj) 
l-Es{u;)Sn{uj)' 



(7) 



where, Ed(uj), Ei^(uj), and Eg {to) are complex error coef- 
ficients, representing the directivity, the reflection track- 
ing, and the source mismatching, respectively. A stan- 
dard method to determine these unknown error coeffi- 
cients is to measure three known reference samples as a 
function of frequency at each temperature. Typically, we 
used a gold film as a short standard, a NiCr film as a load 
standard, and a teflon sheet as an open standard, follow- 
ing the work by Stutzman et al^ By performing this 
standard calibration procedure, both ai{Lo) and '72{lo) 
for all LSCO films were successfully obtained in the fre- 
quency range at least between 0.1 GHz and 10 GHz, as 
shown in Fig.[2l For all films, we confirmed that the skin 
depth 5 at 10 GHz, which was given by cr(a;), was suf- 
ficiently longer than the film thickness t over the whole 
measured temperature range including the vicinity of Tc, 
and that the frequency dependence of both (Ti(w) and 
0-2 (w) well above Tc could be regarded as those in the 
Hagen-Rubens limit of the Drude conductivityii^^ 
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C. Calibration procedure to obtain complex 
fluctuation conductivity 

The standard calibration procedure using three known 
standards is based on the high reproducibihty between 
four measurements including three known reference sam- 
ples and an unknown LSCO sample. In our measurement 
apparatus, the measured Su was reproduced typically 
within 0.05 dB in the magnitude and 0.2° in the phase 
of Sii at 1 GHz. As shown in Fig. 2, this reproducibility 
was enough for us to obtain a{u!) between 0.1 GHz and 
10 GHz for LSCO films far above T^. However, as was 
already reported by our previous paper, we found that 
the present reproducibility was insufficient to obtain (t{uj) 
in the vicinity of Tc- This is because a small unexpected 
difference between the phase of 5'^^'°''* and that of 5'}°'''^ 
gave rise to a large error in <7{lu) in the vicinity of Tc, as 
suggested by Eq. ([6]). 

In order to overcome this difficulty, we modified the 
calibration procedure as follows. First of all, we regarded 
both the magnitude and phase parts of Sii{uj) for the 
LSCO sample, which were measured at a temperature 
(T = To) well above Te, as those of S\°:^'"^, based on the ob- 
servation that the normal-state conductivity of the LSCO 
films was regarded as that in the Hagen-Rubens limit of 
the Drude conductivity. Next, we equated the phase of 
Sfl°'^ to that of at T = To, based on the assump- 

tion that a small difference between them was purely due 
to the experimental error. Finally, we assumed that the 
phase parts of both and 5*}°^"^ were independent of 

temperature in a narrow range between Tc and To. Thus, 
the phase of 5*12°"^* agreed with that of 5'}°'''* at any tem- 
perature from Tc to Tq, as expected for the ideal load 
and short standards. We also confirmed the validity of 
this procedure by the same measurements of conventional 
NbN superconducting thin films as a reference*^ 

If cr(u;) was obtained at each temperature in the vicin- 
ity of Tc through the modified calibration procedure, the 
complex fluctuation conductivity, CTfi(a;,T), is given by 
the following equation, 

<jfi{uj,T) = <j{Lj,T)-a„{T). (8) 

Here, an(T) is the normal-state component included in 
a{u!,T), which agrees with the normal-state conductiv- 
ity far above Tc. Although the change of cr„(T) with 
decreasing temperature was estimated by the linear ex- 
trapolation of dc conductivity from higher temperatures, 
we found that the dynamic scaling behavior of the ex- 
tracted crfl(a;) in the vicinity of Tc was insensitive to the 
temperature dependence of (t„(T). 

IV. RESULTS 
A. Doping dependence in zero magnetic field 

Figure [3] shows the typical frequency dependence of 
a{u!) for the optimally doped LSCO (a;=0.16) and the 



overdoped LSCO (a;=0.20) at several temperatures near 
Tc in zero magnetic field. We observed that both ai{Lu) 
and (72 (t^) in the low frequency limit diverged rapidly 
as the temperature approached Tc from above, suggest- 
ing that the contribution of the superconducting fluctu- 
ations to (j(uj) was evident with decreasing temperature. 
This behavior was more clearly demonstrated by plot- 
ting the phase and the magnitude of afi{Lu){= a{uj) — a„) 
as a function of frequency. Figures lUJa) and (c) shows 
that the plots of the dimensionless phase angle, of tTfl 
versus frequency move to a lower frequency region with 
decreasing temperature, suggesting that a characteristic 
time scale of the superconducting fluctuations is critically 
slowing down. On the other hand, as shown in Figs.|Hb) 
and (d), the magnitude, |cr|, of an increases rapidly with 
decreasing temperature, suggesting a critical divergence 
of the superconducting fluctuations. 

Following the scaling procedures described in Sec. HIl 
we tested whether these behaviors were truly attributed 
to the critical fluctuations of the superconducting order 
or not. Figure [5] shows the results of the data scaling of 
ipcr (upper panels) and |(Tfl| (lower panels) for the LSCO 
films with five Sr concentrations from a:=0.15 to a:=0.20. 
As clearly indicated in Fig. [5l we confirmed that both of 
and |(t| successfully collapsed on to single curves for all 
Sr concentrations, strongly suggesting that the dynamic 
scaling hypothesis was successfully satisfied so that we 
could investigate the change of the critical behavior as a 
function of the Sr concentration along this line. 

The solid and dashed lines depicted in Fig. [5] were 
the 2D and 3D Gaussian scaling functions, respectively, 
which were calculated by Schmidt.'*^ They play roles of 
the guide to eyes for the dimensionality in the critical 
charge dynamics, suggested by Eqs. (4) and (5). We 
found that the optimally doped LSCO (a;=0.15, 0.16) 
showed the 3D-like critical behavior while the overdoped 
LSCO (a;=0.18--0.20) showed the 2D-hkc critical behav- 
ior.'*'* This suggestion was quantitatively confirmed by 
the plots of the scaling parameters, u/q and ao, versus 
the reduced temperature, T/Tc — 1, as shown in Fig. [6l 
Here, we used the values of T^'^^''^ as Tc. In the dynamic 
scaling theory by FFHfi^ the plots of u;ocro(oc are 
expected to be independent of T/Tc — 1 for d—2, while 
they are expected to increase monotonously with T/Tc— 1 
for d=3. Thus, the plots of cjocto in Figs. E^a) and (d) 
clearly indicate that the dimensional crossover from 3D 
to 2D occurred near the boundary between the optimally 
doped region and the overdoped region. 

We also confirmed that luq and cto for the optimally 
doped LSCO with x=0.15 and x=0.16 agreed with 
those for the relaxational 3D-XY universality class with 
1^=0.67 and 2R:!2,'*^ as shown in Figs. Eljb) and (c). On 
the other hand, u/q and ctq for the overdoped LSCO with 
x=0.18-0.20 were found not to show an exponential sin- 
gularity expected for the 2D-XY universality class, but 
rather to show a divergence with an exponent of vz^l.b, 
as shown in Figs. [6lje) and (f). This shows a sharp con- 
trast to the behavior for the underdoped LSCO with 
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a;=0.07-0.14, as already reported in the previous paper3 
We note that the critical behavior with i/z~1.5 ob- 
served for the overdoped LSCO is different from both the 
BKT transition classified in the 2T)-XY universality class 
and the 2D Gaussian fluctuation whose critical exponents 
are v—0.5 and z=2. Thus, we refer this unusual behav- 
ior to the 2D-"U", as an unknown universality class. In 
the next section, we discuss the origin of this unusual 
behavior, in terms of various possible candidates. 

Figure [7] shows a summary of our results about the 
critical charge dynamics of LSCO as a function of the 
Sr concentration. The temperature region where the 
suggested universality class became dominant was deter- 
mined by fitting of the temperature dependence of Wq to 
corresponding models. Together with our previous study 
for the underdoped LSCOj^- we found that two kinds 
of the sudden dimensional crossovers took place across 
the phase diagram. One is the dimensional crossover 
from 2D-XY to 3D-XY between x=O.U and x=0.15, 
while the other is the dimensional crossover from 3D- 
XY to 2D-"U" between x=0.16 and a;=0.18. The first 
dimensional crossover was already reported in our previ- 
ous study^ while the second was observed by this work 
for the first time. 

In the following two subsections, we first investigate 
the critical behavior in the overdoped LSCO in more de- 
tail, in terms of finite size effects. 



B. Finite size effect 

In this study, we used relatively thin LSCO films in the 
overdoped regime, as listed in Table I. This is because, 
for the reliable data of (Tfl(a;) to be obtained by the mi- 
crowave broadband technique, the sheet resistance just 
above Tc is needed to be larger than about 2 fi. This 
condition requires us to use thinner films for less resis- 
tive overdoped LSCO. However, the use of thinner films 
leads to a suspicion that the observed 2D-like critical 
behavior for the overdoped LSCO is caused by a finite 
film thickness. Indeed, we have demonstrated that the 
dynamic fluctuations in the superconducting NbN thin 
films showed the 2D Gaussian fluctuations, which were 
introduced by the divergence of ^ beyond the film thick- 
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ness.— 

In order to investigate this problem, we measured the 
films with different thickness near the Sr concentration 
where the second dimensional crossover occurred. Fig- 
ures El and [9] show the scaled data of ^o- for a:=0.17 
(t=115, 140 nm) and x=0.18 {t=60, 240 nm), respec- 
tively. 

First of all, as shown in Fig. [H the behavior of 0(t(w) 
scaled by coq for a;=0.17 showed the 2D-like behavior for 
the thicker film (t=140 nm), and the 3D-like behavior 
for the thinner film (t=115 nm). It is clear that this 
change of the dimensionality cannot be explained by the 
finite size effect, in contrast to the previous results on 
the NbN films. Rather, this result implied a possibility 



that the intrinsic dimensional crossover occurs suddenly 
near x=0.17. 

The behaviors of 00- for a;=0.18 agreed with the 2D- 
like critical charge dynamics even for the thicker film 
{t—240 nm), as shown in Fig. [5) By comparing this 
result with that of the film with a;=:0.16 (t=140 nm) 
which showed the 3T)-XY critical charge dynamics, it 
provided us the following important insight. As shown 
in Fig. [Tni the 2D-like critical fluctuations for the thicker 
film could be observed in the temperature range from 
1.006Tc to 1.02rc. If we consider that the 2D behavior 
is due to the finite-thickness effect on the 3D-XY criti- 
cal behavior, ^ should be larger than the film thickness 
(t=240 nm) below 1.02Tc. This requires that is larger 
than 20 nm. On the other hand, the 3D-XY fluctua- 
tions for x=0.16 was observed in the temperature range 
from 1.004Tc to 1.02Tc. This means that ^(r) is not 
larger than t=140 nm above T=1.004Tc, requiring that 
^0 is less than 3.5 nm. Thus, a difference of between 
a;=0.16 and a;=0.18 need to be larger than 6 times in 
magnitude. This seems to be too large to be explained 
by the change of Tc or the superconducting gap with 
hole doping. Therefore, the observed behaviors in (f>„ for 
these films with different thickness never supported that 
the 2D-like critical behavior near a;=0.18 was attributed 
to the finite thickness effect. The same conclusion was 
also obtained from the scaled data of \a\. 



C. Effects of external magnetic field 

Another way to investigate a finite size effect in the 
critical charge dynamics is the use of an external mag- 
netic field parallel to the film. If a magnetic field larger 
than the lower critical field was applied to a type-II su- 
perconductor, vortices can penetrate into the supercon- 
ductor. Since the vortex pinning effect is negligibly small 
near Tc, vortices are distributed uniformly over the sam- 
ple. At the core of each vortex, the growth of the super- 
conducting order parameter is strongly supressed, sug- 
gesting the divergence of ^ due to the superconducting 
fiuctuation is cut off by a finite distance between vortices. 
Thus, the external magnetic field can be used as a tun- 
able probe to study the finite size effect, since the mean 
distance between vortices, £v, is proportional to ^J^q/B. 

Of course, it is crucially important to distinguish the 
critical charge dynamics under a magnetic field from the 
vortex dynamics in the mixed state. For this purpose, we 
measured four LSCO films (a;=0.07, 0.16, 0.17 and 0.18) 
by applying magnetic field perpendicular to the Cu02 
planes up to 1 T systematically and investigated the ob- 
tained results very carefully. 

Figure [TT] shows the scaled data of 4)cj and for the 
underdoped LSCO (x=0.07) at several constant magnetic 
fields. We found that both <j)a and \a\ were scaled suc- 
cessfully up to i?=0.1 T, while neither of them were col- 
lapsed on to single curves at B=l T. This suggests that 
the contribution of the vortex dynamics is almost negli- 
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gible up to at least B=0.1 T, whereas it cannot be ne- 
glected at B=l T. As shown in Figs. [TWa) and (b), we 
confirmed again that the behavior of loq and tro obtained 
at B—0 T excellently agreed with the BKT theory in 
the temperature range from 1.02Tc to 1.2Tc. Note that 
ujq^ cx ctq oc ^^rp in the relaxational 2D-XY universality 
class. Here, the BKT correlation length, ^kt, shows an 
exponential divergence in a critical region, as follows. 



^KT^Coexp 6/VT/r, ~ 1 



(9) 



where, & is a numerical constant, which was determined 
to be ~ 0.2 in this study. 

The plots of the correlation length, ^, given by 



as a function of 1/ y^T/Tc — 1 were found to be very use- 
ful to understand the finite size effect due to the external 
magnetic field. As is shown in Fig. [TWc). the divergence 
of ^ in the vicinity of Tc was suppressed with increas- 
ing magnetic field. This behavior can be understood by 
considering two kinds of the characteristic length scale; 
the BKT correlation length, ^kt, and the mean distance 
between vortices, £v In the temperature region not very 
close to Tc, where ^kt ^ ^v, the logarithmic interaction 
between a vortex and an anti-vortex, which is essential 
to the BKT transition, is not screened by the externally 
induced vortices. Thus, the exponential divergence due 
to ^KT is observed. On the other hand, in the vicinity 
of Tc, where ^kt > ^u, the externally induced vortices 
can screen out the interaction between vortex and anti- 
vortex, leading to a crossover of the characteristic length 
scale from ^kt to £y. Thus, the growth of ^ is suppressed 
as f approaches iy. We confirmed that the saturated 
value of 1^ in the limit of Tc did not exceed £y at each 
magnetic field, by using a reasonable assumption that 
^0^3 nm. In summary, the result of Fig. [T^is also con- 
sistent with the BKT picture for the fluctuation of the 
x=0.07 sample. 

Figures [T3] and [U show the results of dynamic scal- 
ing analyses for the optimally doped LSCO (a:=0.16) at 
several constant magnetic fields. We found that the data 
collapse of the scaled (l)a and \a\to single curves were suc- 
cessful up to B=l T, suggesting that the contribution of 
vortex dynamics was almost negligible in the measured 
magnetic field region. In addition, all the results clearly 
indicated that the 3D-XY critical charge dynamics was 
maintained up to B—1 T. Thus, the external magnetic 
field seemed not to give any change to the critical charge 
dynamics near Tc, in contrast to the case of the BKT 
transition. This is because the power-law like divergence 
of £,{T) in the 3D-XY critical charge dynamics is much 
slower than the exponential divergence of ^kt(T). In- 
deed, by assuming that ^o~3 nm and ^—0.67, we could 
estimate that f(T) for this film (a;=0.16) was less than 
31 nm in the observed temperature range at B—1 T. 
Since £y is estimated to be ~ 45 nm at B=l T, we can 
safely conclude that ^(T) is always smaller than £y even 
at B=l T, strongly suggesting that the critical behavior 
is never affected by the external magnetic field at least 
up to B=l T. 



Figures [15] and [16] show the results of dynamic scaling 
analyses at several constant magnetic fields for the over- 
doped LSCO film (x=0.18 and t=60 nm). As was already 
pointed out, the critical charge dynamics of this film in 
zero magnetic field was represented by an unknown uni- 
versality class in the 2D system (2D-"U"). Figure [T5I 
clearly indicates that the scaled data of both ipa- and \<t\ 
up to B=l T were very similar to those at B=0 T. In 
addition, the temperature dependence of the scaling pa- 
rameters at all magnetic fields except for B=l T excel- 
lently agreed with each other, as shown in Fig.]!^] These 
results strongly suggest that the small applied magnetic 
field (at least up to 0.3 T) did not affect the critical charge 
dynamics of the overdoped LSCO. 

Thus, for both a;=0.16 and 2;=0.18, the magnetic field 
dependence of the critical temperature, Tc{B), which was 
experimentally determined through the dynamic scaling 
analysis at finite magnetic field below _B=0.3 T, can be 
regarded as the same as the temperature dependence of 
the upper critical field, i?c2(T), as shown in the inset of 
Fig. [T7J Note that Bc2 is not a phase transition point 
in practice. There is only a first-order transition of vor- 
tices at a lower field. However, the concept of Bc2 as 
the mean-field critical field still makes sense. This pro- 
vides us another route to estimate v and ^o, by using the 
following relationship, 



Sc2(T) 



27rC(T)2 27r 



(1 - T/Tco)^ 



(10) 



where Tco is the critical temperature at B~0 T. As shown 
in the main panel of Fig. 1171 we confirmed that the 
plots of the applied magnetic field, B, as a. function of 
(l-Tc(5)/Tco) for the optimally doped LSCO (a;=0.16), 
agreed with a line of Bc2(T) = (/)o/2<(T)2 with i^=0.67 
and ^0^3 nm, as suggested by the dynamic scaling anal- 
yses for a;=0.16. This supports the idea that the plots of 
_B as a function of (1 — Tc{B)/Tco) can be regarded as 
those of Bc2(T) as a function of (1 - T/Tco). 

Based on this idea, we estimated that and 
^o'~-"3 nm for the overdoped LSCO (a;=0.18). As shown in 
Fig. [TT] it is clear that the value of v for a;=0.18 is larger 
than the values for the GL theory (i/GL=0.5) and the 
3D-XF universality class (i'3dxy=0.67). We found that 
^(T) increased up to 70 nm in the measured temperature 
range at B—0.3 T, by using the estimated values of v and 
^0- Since £y was estimated as 83 nm at 5=0.3 T, ^ was 
always smaller than £y within the measured temperature 
range at least up to i?=0.3 T, which confirmed that the 
critical charge dynamics, not the vortex dynamics, was 
observed. Together with the value of i/z{^ 1.5) obtained 
at B=0 T, we obtained that 2:^1.7 for a;=0.18, which was 
slightly smaller than the well-known value for the relax- 
ational dynamics (z«2). Thus, the experimental results 
for a;=0.18 suggested that the 2D-"U" universality class 
had the critical exponents of j/~0.9 and 2;~1.7. 

In contrast to the results for a::=0.16 and 0.18, the 
behavior of (j)a and \a\ for a;=0.17 (the thicker film, 
t=140 nm) was largely dependent on the applied mag- 
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netic field. Up to 0.1 T, the scaled data of 4>a and |(t|, 
and the temperature dependences of the scaling parame- 
ters were similar to those at i?=0 T, as shown in Figs.fTSl 
and I19[ respectively. By assuming that was about 
3 nm, which was the same value as those used for x—Q.\Q 
and 0.18, ^(T) at 1.02rc was estimated to be 100 nm 
for the 2D-"U" critical charge dynamics. Since ly was 
about 140 nm at i?=0.1 T, ^(T) did not exceed up 
to _B=0.1 T. Thus we confirmed that the contribution 
of vortex dynamics was negligible and that the 2D-"U" 
critical charge dynamics retained up to 0.1 T. However, 
above 0.3 T, we found that the values of (p^ioj /ujq) in 
the high frequency limit clearly decreased with increasing 
magnetic field, and also that the frequency dependence of 
|cr(w/a;o)| was weakened, as shown in Fig. [TBI Within the 
dynamic scaling theory, ^^'^^ these results seem to suggest 
that the dimensionality in the critical charge dynamics 
changes from 2D to 3D by applying magnetic field. Fig- 
ure [19] also shows that the temperature dependence of 
ctq Wo (oc ^^~'^) was changed markedly to the 3D-hke above 
0.3 T. However, we could not clearly confirm the 3D-XF 
critical charge dynamics, since ujq and ctq could be fitted 
by the critical exponents of the 3D-Xy universality class 
only in a very narrow region between l.OlTc and 1.03rc. 
This suggests a possibility that another contribution af- 
fects the critical behavior above 0.3 T. This possibility 
was also suggested by the fact that the scaled data of 
(pa at B—1 T were more scattered than those at _B=0 T, 
leading to poorer fitting results of loq and ctq to the 3D- 
XY model. As one possible candidate, we suspect that 
the vortex dynamics dominates the data. Indeed, the 
estimated ^(r)~100 nm for the 2D-"U" critical charge 
dynamics at 1.02rc exceeds £t,~80 nm at 0.3 T. How- 
ever, we found that the scaled data of (j),, for x=0.17 
above 0.3 T was qualitatively different from those ob- 
served for a;=0.07 at 1 T, as shown in Figs. [TT] and [TH] 
In addition, if we use the 3D-XY critical exponent, ^(T) 
in the measurement temperature range is estimated at 
most ^--^60 nm, which is smaller than €t,~80 nm at 0.3 T. 
Therefore, there is another possibility that the critical 
charge dynamics changed from the 2D-"U" to the 3D- 
XY universality class, even though it remains unclear at 
present why the critical charge dynamics for a:=0.17 was 
largely dependent on the applied magnetic field. 

In summary, the systematic studies of the effects of 
the external magnetic field on the critical charge dynam- 
ics confirmed that our determination of the universality 
class, which was classified into 2T)-XY for the under- 
doped region, iD-XY for the optimally doped region, 
and 2D-"U" for the overdoped region, was plausible. 



V. DISCUSSION 

As was already pointed out, we discovered three differ- 
ent universality classes in the phase diagram of LSCO, 
representing the existence of two kinds of crossover lines 
across the phase diagram. One is the dimensional 



crossover from 2T)-XY to ZD-XY between a;=0.14 and 
x=0.15, while the other is the crossover from ZT)-XY 
to 2D-"U" between a;=0.16 and x^O.lS. The study of 
the effects of thickness and magnetic field dependence 
supported that these features were not caused by the 
finite size effects, but were intrinsic to the superconduct- 
ing phase transition in LSCO. In this section, we discuss 
possible origins of these changes of the universality class 
with the Sr concentration from various points of views. 



A. The 2Y)-XY critical charge dynamics in the 
underdoped region 

First, we discuss the 2D-XY critical charge dynamics 
in the underdoped region {x = 0.07 ^ 0.14). This pos- 
sibility was already pointed out by our previous studyi^ 
In this work, we confirmed that the exponential diver- 
gence of CKT(r) for x=0.07 was suppressed by applying 
small magnetic field, in agreement with the BKT the- 
ory which considered the effects of externally induced 
vortices4^ These new results strongly support that the 
critical charge dynamics in the underdoped region was 
classified into the 2D-XY universality class. What is 
most significant is that the phase fluctuations of the su- 
perconductivity order parameter, which could be directly 
determined in our methods^^ survived only up to much 
lower temperature (at most ^1.4Tc) than the closing tem- 
perature of the pseudogap. As shown in Fig[71 this shows 
a sharp contrast to the Nernst experiments by Ong and 
Wang4^ Our results indicated that the anomalous pseu- 
dogap phenomena cannot be explained only by the classi- 
cal phase fluctuations of superconducting orders, in con- 
trast to the early prediction by Emery and Kivelson.'^ 
A similar conclusion was also obtained by a theoretical 
consideration on "cheap" vortices where the vortex core 
energy is very small (~ k^Tc) and both of the amplitude 
and phase fluctuations are controlled by the same energy 
scale, ksTci^ 

In the layered superconductor, it is suggested that the 
finite interlayer coupling effect enhances the true critical 
temperature, Tc, as follows. 



BKT 



= 1 



Inyi/A 



(11) 



where Tbkt, A are the BKT transition temperature 
and the anisotropy ratio of the interlayer coupling to 
the intralayer coupling, respectively i^°i^^ According to 
Eq. (fTT|) . the 2T)-XY critical behavior is no longer ob- 
served in the vicinity of Tc and that the system should 
exhibit the 3D critical behavior, except for the case that 
the interlayer coupling is negligibly small (A 0). How- 
ever, it should be pointed out that Eq. (fTTjl overesti- 
mated Tc- In fact, Eq. pT|) was derived by considering 
the critical radius, Vc, where the energy of the indepen- 
dent vortex pair balances with that of the vortex ring 
(See Fig. [20]). In Refs. HO and IsH. Vc was estimated by 
neglecting the contribution of the vortex core energy, Ec, 



9 



compared with the logarithmic interaction term. How- 
ever, this is incorrect, which was aheady pointed out in 
the original paper of the BKT theory,— and also in a re- 
cent theoretical review paper more clearly42. Particularly, 
according to the latter paper, Ec is expected to be at least 
comparable to the logarithmic interaction term even for 
"cheap" vortices. Thus, it is clear that Eq. (fTT|) overes- 
timates the enhancement of Tc owing to the interlayer 
coupling effect. In order to discuss the possibility of the 
vortex ring excitation in the layered cuprate supercon- 
ductor more precisely, we need the accurate estimation 
of Ec of "cheap" vortices, which remains an open issue. 

Even if we consider a somewhat unrealistic situation 
where Ec is negligible, our results merely suggest that 
the interlayer Josephson coupling is not yet developed at 
least above the minimum temperature where the critical 
behavior was successfully observed (~ l.OlTc for x=0.Q7, 
as shown in Fig. [7]). Indeed, we found that the effective 
thickness of the superconducting sheet was comparable 
to the distance between the Cu02 planes, suggesting that 
the fluctuating superconductivity appeared in each Cu02 
plane which remained nearly decoupled within the mea- 
sured temperature range.— We emphasize that our re- 
sults do not rule out the development of the interlayer 
Josephson coupling, which will finally lead to the 3D crit- 
ical behavior. The Josephson coupling is expected to be 
developed in the further vicinity of Tc, which we could 
not approach experimentally. Thus, our observation of 
the 2D critical behavior does not contradicts the general 
property of a layered superconductor. 



B. Two crossovers in the phase diagram 

The most important findings are the two crossovers, 
that is, the crossover from 2D-XY to 3D-XY near 
x=0.14^ and the other crossover from 3D-XY to 2D- 
"U" near a;=0.16 in the phase diagram. First, we focus 
on the first crossover. The sudden change of the univer- 
sality class from the 2D-XY to 3D-XY around at x—O.lA 
is discussed in terms of several theoretical models. 

In the U{1) mean field theory of the t — J model^i^ 
Tc is controlled by the Bose condensation of holon in the 
underdoped region, and the fermion pairing in the over- 
doped region. This means that a qualitative change of 
the condensation nature is expected to take place at the 
optimal doping. Our results were consistent with this 
picture in the sense that the critical charge dynamics 
changes near the optimally doped region. However, it 
seems to be difficult to explain the change of the dimen- 
sionality. In the theory, the weak 3D interlayer hopping 
is assumed to obtain a finite transition temperature of 
the Bose condensation of holons. As was suggested by 
Lee and co-workers?^ the BKT transition for the fermion 
pairing and the boson condensation can occur simultane- 
ously where the gauge field becomes massive due to the 
Higgs mechanism. Thus, the phase transition lines for the 
fermion pairing and the boson condensation, which were 



predicted in the C/(l) mean filed theory of the i— J model, 
become the crossover lines and only the superconducting 
transition remains to be the real BKT transition. This 
suggests that we should observe only the 2D-XY univer- 
sality class in the whole range of the phase diagram. In 
other words, the anomalous change of the dimensionality 
in the critical behavior with carrier doping cannot be ex- 
plained only by the early framework in the t — J model 
suggesting the fermion pairing and the boson condensa- 
tion. The so-called confinement-deconfinement transition 
due to the strong coupling gauge field^^'^^ is another can- 
didate for the change of the critical behavior. Although 
a confinement-deconfinement transition does not always 
imply the change of the dimensionality, it is possible that 
a 3D Fermi liquid phase is realized as a consequence of 
the confinement of fermions and bosons, in which the 
3D-XY critical behavior arises. 

On the other hand, in the hidden QCP scenario, this 
sudden change of the universality class can be understood 
as a classical-quantum crossover which may occur with x 
approaching the hidden QCP located beneath the super- 
conducting dome. A spatial and time correlation of the 
quantum fluctuation of the hidden order is strong near 
the QCP. Therefore, provided that the quantum fluctu- 
ation couples to the thermal fluctuation of the super- 
conducting order, the critical dynamics around the su- 
perconduting transition can be affected by the quantum- 
fluctuation effects in the vicinity of the QCP. In this pic- 
ture, it is also predicted that the quantum fluctuation 
effects should become less prominent as x increases away 
from the hidden QCP. From this point of view, the be- 
havior in the overdoped region reported in this paper is 
very crucial to examine these theoretical models. 

In the overdoped region, we found the second crossover 
from -SD-XY to 2D-"U" near a;=0.16, in addition to the 
first crossover near a;=0.14. This is highly unexpected, 
since the general trend is that the transport properties of 
the high-Tc cuprates become less anisotropic with carrier 
doping. — The t — J model based interpretation encoun- 
ters a difficulty, since it does not expect any further di- 
mensional crossover to take place. Even if a 3D Fermi liq- 
uid phase was realized by a confinement of fermions and 
bosons, there is no reason why the 3D phase is realized 
only in a narrow range near the optimally doped region. 
On the other hand, the theoretical models assuming the 
hidden QCP can naturally explain our finding that the 
critical charge dynamics at least down to T^l.QlTc was 
essentially 2D all over the phase diagram except for a 
narrow range near the optimally doped region. This re- 
sult was in correspondence with the expectation that the 
quantum fluctuation effects are prominent only in the 
vicinity of the QCP. 

More precisely, the quantum fluctuation effects at finite 
temperatures are prominent only in the quantum critical 
region, which will be described below. In general, in the 
2D-QPT scenario, the hidden order of the QPT becomes 
long-range only at zero temperature.^* For instance, sup- 
pose the Heisenberg antiferromagnet with the 2D square 
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lattice, which corresponds to an undoped Cu02 planci^ 
The magnetic long-range order is present only for g < gc 
at T = 0, where g and gc are a key parameter to con- 
trol a QPT and a critical parameter corresponding to the 
QCP, respectively. At T = 0, the correlation length of 
QPT diverges as ^ \g — gd '^ in the vicinity of gc, where 
is a static quantum critical exponent. On the other 
hand, at finite temperatures, there are three characteris- 
tic regions, each of which is separated by two crossover 
lines, T — ±{g — gcY^ , where z is a dynamic quantum 
critical exponent, as shown in Fig. 1211 Here, we used 
the same description for each region as that proposed by 
Sachdevfi^ since the commonly used description such as 
"quantum disordered" seemed to be somewhat mislead- 
ing. 

In the low T region on the magnetically ordered side 
[0 < r < {gc — gY^]i both of and the correlation time, 
rg, decrease exponentially as T increases, suggesting that 
the quantum fluctuation effect is extremely small (line(a) 
in Fig. [2T|). Thus, in this region, the system is regarded 
as classical at finite temperatures. In the quantum criti- 
cal region [T > {gc- gY"" and T > {g-gcY^}^ and tq 
decrease as ~ 1/T (line(b) in Fig. [2T|) . This means that 
large quantum fluctuations survives even at high temper- 
ature which might appear in various properties, as well 
as thermal fluctuations. Finally, in the low T region on 
the quantum paramagnetic side [0 < T < {g — gcY^']i 
^g exponentially reaches a constant value as T — > 
(line(c) in Fig.[2T|). The value of ^g at T = behaves as 
Iff ~ ffcT^- This region can also be regarded as a "spin 
liquid" state,— since there is no long-range order even 
at T = 0. In this region, although quantum fluctuations 
surviving at T = contribute to the oscillatory spin cor- 
relation, the long-range time correlation is exponentially 
relaxed by thermal fluctuations. Thus, the behavior of 
Tg at T > is rather dominated by the classical relax- 
ation. In summary, the schematic behavior of ^g for the 
2D Heisenberg antiferromagnet is shown in Fig. [21] 

When this concept is applied to discuss the phase di- 
agram of the cuprates, the parameter, g, is related to 
the carrier concentration, n. As holes are doped into 
the Cu02 plane at T = 0, it is expected that the long- 
range order of the spin correlation is broken at a critical 
value of the hole doping, ric, which is corresponding to 
gc- Thus, one can consider that 5 is a positive function 
of the hole doping and that g{nc) = gc- In the quantum 
critical region around Uc, the antiferromagnetic quantum 
fluctuation is considered to decay slowly as ^ l/T. It is 
already well-established that the antiferromagnetic quan- 
tum fluctuation favors the d-wave superconductivity4i^ 
This suggests the possibility that a spatial and time cor- 
relation of the superconducting order can be enhanced 
by the quantum fluctuation effects. Provided that the 
interlayer correlation is also enhanced, the observed di- 
mensional crossover from 2D to 3D can be understood as 
a consequence of this quantum-fluctuation effects. Note 
that, this effect is limited to the quantum critical re- 
gion. In other two regions, rg for the quantum fluctu- 



ation decays exponentially with increasing temperature, 
suggesting a purely classical dynamics. Thus, the exper- 
imental result that the 3D behavior was observed only in 
the vicinity of the optimally doped region showed good 
qualitative agreement with this QCP scenario, as shown 
in Fig. da 

Note that the 2D antiferromagnetic order discussed 
here is just one of the candidates for the possible QPT. 
The scenario is essentially based on a hidden order in the 
2D system. The 2D order is present only at T = 0, and 
becomes short-range at finite temperatures, independent 
of g- Therefore, no experiment at finite temperatures 
can determine Uc directly. This means that the disap- 
pearance of the 3D Neel order by hole-doping does not 
always imply the position of the QCP in the 2D antifer- 
romagnet. Rather, a scaling analysis of the wave vector- 
and frequency-dependent magnetic fluctuations seems to 
be useful to search the position of the QCP. Indeed, the 
neutron scattering study using an LSCO single crystal 
suggested that the QCP might be hidden near x—O-U^ 
which agrees well with our results. 

At present, this scenario seems to be compatible with 
our results most successfully. The important requirement 
is that the quantum fluctuation around the hidden QCP 
couples to the thermal fluctuations of the superconduct- 
ing order near the superconducting transition. Thus, a 
crucial question is what is the most plausible candidate 
for a hidden order among various proposals 
A very recent discovery,— of the 4ao wide unidirec- 
tional electronic domain which was detected by using a 
tunneling-asymmetry imaging technique in STM, may be 
also a candidate for such hidden orders. Unfortunately, 
we could not find a clear answer to this question in this 
work. In addition, there is a more complicated problem 
in the neighborhood of the Mott insulator-to-metal tran- 
sition. It is predicted that the coupling between the order 
parameter of the QPT and the low-energy fermions pro- 
vides a more complex behavior with two dynamic quan- 
tum critical exponentsi^i^'^ In LSCO single crystals, 
the insulator-to-metal crossover was suggested to occur 
near x=0.16, by the in-plane resistivity measurements 
down to T=0.65 K using a 61 T pulse magnetic fieldi^ 
An open issue is whether this behavior is attributed to 
the Mott transition or not. 



C. Possible origin on the 2D-"U" critical charge 
dynamics 

The critical charge dynamics observed in the over- 
doped region was found to be different from the well- 
known universality classes {2'D-XY, 3D-XY, or the 
Gaussian fiuctuation). Here, we discuss this 2D-"U" crit- 
ical behavior in details. 

First of all, we consider the effect of disorder due to the 
Sr substitution, which is expected to be more prominent 
with increasing Sr concentration. As was emphasized re- 
peatedly, the successfully scaled data confirmed that the 
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observed 2D- "U" critical behavior could not be explained 
by the distribution of due to the disorder effect, be- 
cause our method to obtain the data collapse of 0o- and 
\a\ is sensitive to the breakdown of the dynamic scaling 
hypothesis, which can be caused by a broad distribution 
of Tc due to disorders. Figure [23] shows an example of 
the breakdown of the dynamic scaling due to the distri- 
bution of Tc- This is 0o-(^) scaled by luq of the thicker 
film (t=120 nm) with x = 0.18. As shown in the inset of 
Fig. [221 dc resistivity showed a broader transition width 
than that of the thinner film {t—60 nm), which yielded 
the successfully scaled cta (w) . The data of the thicker film 
does not show the dynamic scaling behavior, suggesting 
a broad distribution of Tc due to disorders results in the 
breakdown of the dynamic scaling hypothesis. 

In addition to such an experimental justification, the 
so-called Harris's criterion on the critical phenomenon 
also tells us that the weak disorder is irrelevant to the 
2D-XY or SB-XY critical behavior.^. On the other 
hand, the strong disorder can change the critical behav- 
ior to a different one. However, even for the strong disor- 
der, it seems difficult to change the dimensionality in the 
critical charge dynamics. In other words, the observed 
dimensional crossover from 3D to 2D is a very robust 
one, which cannot be explained by the disorder effect. 

Meanwhile, the strong disorder effect is one of the pos- 
sible candidates for the origin of the anomalous critical 
exponents. The disorder effect on the 2D system has been 
investigated by the 2D random gauge XY model, which 
introduced the disorder effect as the random variables in 
magnetic bond angles uniformly distributed in the range 
from —rir to rir^ When the disorder strength, r, is equal 
to 1, the model corresponds to the standard XY gauge 
glass modelr^ while the standard XY model is recovered 
for r = 0. Interestingly, a recent Monte Carlo simulation 
of the 2D random gauge XY model suggested that a non- 
BKT type transition with v ^ l.l occurred at low tem- 
peratures for a stronger disorder strength than a critical 
value of r, while the standard BKT transition occurred at 
low temperatures for a weaker disorder strength;^ The 
calculated critical exponent is not so different from our 
result {i^ ~0.9) for a;=0.18. Thus, this is one possibil- 
ity for the 2D-"U" universality class. However, a better 
agreement of the critical exponent is obtained by consid- 
ering the frustration effect on the 2D system. We will dis- 
cuss this later in terms of the frustrated 2D-XY model. 

Next, in the t — J model assuming the BKT transition 
for the fermion pairing and the boson condensation, it is 
expected that the 2D-XY critical behavior is observed 
in a whole region of the phase diagram. In practice, the 
BKT transition cannot be observed unless the BKT cri- 
terion, A_L ^ L, is satisfied,— where A^(= X^/Dg) and 
L is a screening length and a lateral size of a supercon- 
ducting sheet, respectively, and A and Dg is a bulk pen- 
etration depth and a thickness of the superconducting 
sheet, respectively. In the underdoped region, this crite- 
rion is satisfied, indeed. For overdoped films, A_l was 
estimated to be 10 ^ 100 mm, which is sufficiently larger 



than L 1 mm, showing the criterion is satisfied also 
for overdoped samples. Therefore, the 2D-"U" critical 
behavior was not induced by the failure of the BKT cri- 
terion. Thus, the 2D-"U" feature cannot be explained by 
the t — J model based theories. 

In a model assuming that a Fermi liquid state is real- 
ized in the heavily overdoped region, one may expect that 
a conventional Gaussian fluctuation rather than the crit- 
ical fluctuation becomes more prominent with increasing 
the Sr concentration. From this point of view, a crossover 
from the 2D-XY critical behavior to the Gaussian fluc- 
tuation can be a candidate for the origin of the observed 
anomalous critical exponents. In order to examine this, 
we estimated 1 /r for each film (t is the correlation time 
of the fluctuation), by comparing the scaled data of 4>a-{^) 
with the Gaussian scaling function."'^ In the GL theory 
assuming the Gaussian fluctuation, it is well known that 
the divergence of r near Tc is given as follows. 

" 16kB\T-Tc\' ' 

Surprisingly, we found that all the results of 1/r for 
x=0.07, 0.16 and 0.18 were much smaller than I/tgl by 
two or three orders of magnitude, as shown in Fig. 1241 A 
similar slow critical charge dynamics has been reported 
for the nearly optimally doped YBa2Cu30y thin film.— 
These results show a sharp contrast to the similar plots 
of l/r for NbN thin films, which almost agreed with 
l/'J'GL ( also shown in Fig.[24|). Therefore, the observed 
dynamics is indeed attributed to the critically slowing 
down phenomenon in the critical region, and it cannot 
be explained by a crossover to the Gaussian fluctuation. 
Interestingly, such a slow dynamics is also expected in 
the low T regions both on the magnetically ordered and 
the quantum paramagnetic sides around the QCP, be- 
cause of the opening of an energy gap in the excitation 
spectrum.— 

In terms of the QCP scenario, we believe that the spin 
liquid state in the low T region on the quantum param- 
agnetic side seems to provide key informations to under- 
stand the origin of the 2D-"U" critical behavior. When 
a hole is doped in the Cu02 plane, the spin on a Cu 
site is combined to the spin of the doped hole to form 
the Zhang-Rice (ZR) singletj^ giving rise to the spatial 
modulation of the nearest-neighbor exchange interaction, 
J, to a square lattice in the Cu02 plane. In other words, 
the formation of the ZR singlet by doped holes introduces 
the frustration between the spins in the 2D antiferromag- 
net, as pointed out by Anderson.^- In particular, mobile 
ZR singlets can give rise to the dynamic spin frustration. 
This dynamic or static frustration effect has been consid- 
ered to play an essential role in stabilizing the spin liquid 
or glass state.— 1^ Such a frustrated state might be re- 
lated to the 2D-"U" critical behavior in the overdoped 
region. Thus, we need to consider more explicitly the in- 
fluence of the frustration effect to the superconductivity 
realized in the quantum paramagnetic side. 

The static frustration effect in the XY spin systems 
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has been extensively investigated by the so-called frus- 
trated XY model^ both theoretically^ and experimen- 
tally.— Interestingly, the study of the critical behavior in 
frustrated XY spins seems to suggest the possibility of 
a new universality class, which is clearly different from 
the 2D-XY universality class,- - Indeed, many numerical 
simulations for the fully frustrated XY model suggested 
that 1^=0.8 - 1.0,®^, which shows a good agreement with 
our result (;/~0.9) for x=0.18. Although the explicit cor- 
respondence between the frustrated XY spin system and 
the superconductivity realized in the spin liquid state is 
an open issue to be resolved, the resemblance between 
the critical behaviors of both systems suggests the im- 
portance of frustration effect in the overdoped region. 



VI. CONCLUSION 

We investigated the critical charge dynamics of the su- 
perconducting to the normal-state transition for LSCO 
thin films with a wide range of the Sr concentration 
(a;=0.07 to 0.20), through the dynamic scahng analysis 
of the fluctuation-induced microwave conductivity. We 
discovered that the critical charge dynamics in the phase 
diagram of LSCO was classified into the following three 
different universality classes; (i) the 2D-XY universality 
class for the underdoped region {x—0.07 to 0.14), (ii) the 
3D-XY universality class for the nearly optimally doped 
region (j;=0.15 and 0.16), (iii) the 2D-"U" (z/ ~0.9, z -2) 
universality class for the overdoped region (a;=0.17 to 
0.20). We confirmed that the anomalous 2D critical be- 
havior for the overdoped region is not induced by the 
finite size effect but is intrinsic to the overdoped LSCO. 

Our results indicated that there are two kinds of di- 
mensional crossovers in the phase diagram of LSCO, sug- 
gesting that the dimensionality in the critical charge dy- 
namics changes twice with hole doping. In other words, 
except for the 3D-XY critical charge dynamics observed 
only near the optimally doped concentration, we found 
that the critical charge dynamics of LSCO was essen- 
tially two-dimensional. This suggests that the fluctuating 
high temperature superconductivity initially appears in 
each Cu02 plane, preceding the development of the inter- 
layer Josephson coupling. Our results for the underdoped 
LSCO also indicated that the 2D-XY critical charge dy- 
namics could be observed at most up to ~ 1.4 T^,, which 
was much lower than a closing temperature of the pseu- 
dogap. Thus, some of the early naive pictures about 
the pseudogap and the critical fluctuations in high-Tc 
cuprates need to be modifled. 

We discussed the possible origins of the anomalous di- 
mensional crossovers. It is concluded that neither the 
weakened anisotropic property with hole doping nor the 
conflnement-deconfinement transition discussed in the 
t — J model can explain the second dimensional crossover 
from the 3D-XY to the 2D- "U" , which occurs between 
x=0.16 and 0.17. On the other hand, the so-called hid- 
den "QCP" scenario, where the quantum critical fluctu- 



ation near the hidden QCP couples to the superconduct- 
ing fluctuation, was able to explain the existence of two 
kinds of dimensional crossovers qualitatively. This sce- 
nario also seems to provide a key infromation to under- 
stand the origin of the 2D-"U" critical behavior, whose 
critical exponents were found to be similar to those of the 
frustrated XY model. The importance of the frustration 
effect in the spin liquid state is implied 

At present, the scenario based on a hidden order of 
QPT in the 2D system seems to be the most reasonable 
candidate for the anomalous change of the universality 
class with hole doping. However, the hidden order which 
provides the direct evidence for QCP has not been yet 
discovered. If the hidden order is formed in the 2D sys- 
tem, it will be almost impossible to be found at finite 
temperatures. Thus, the dynamic scaling analysis includ- 
ing the quantum critical fluctuations plays an important 
role to settle this issue. 
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FIG. 2: (a) Frequency dependence of (Ti for LSCO thin films 
with x=0.07 to 0.20 at T^lOO K. (b) Frequency dependence 
of a2 normalized by ai at 0.1 GHz for the same films. For all 
films, (71 is independent of frequency and a2 is much smaller 
than (Ti, suggesting that the charge dynamics in the normal 
state can be safely regarded as in the Hagen-Rubens limit of 
the Drude conductivity. 



FIG. 3: Frequency dependence of (a) ai and (b) a2 of the 
a;=0.16 film at 7=36.0-37.6 K, and (c) ai and (d) C72 for 
the x=0.20 film at r=33.2-36.5 K. All temperatures are just 
above Tc. Both cri and <72 showed a remarkable enhance- 
ment in the low frequency limit as temperature approaches 
Tc. The insets show the temperature dependence of dc resis- 
tivity. Each circle represents the temperature of each data in 
the main panels. 
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FIG. 4: Frequency dependence of (a) the phase, 4>a, and 
(b) the magnitude, \a\, of the fluctuation induced complex 
conductivity for the a;=0.16 film , and (c) 4>a and (d) \a\ for 
the a;=0.20 film. 
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FIG. 5: Scaled data of the phase (upper panels) and the magnitude (lower panes) of afi(w) for the LSCO films with various 
carrier concentrations. The phase and the magnitude is normalized by 7r/2 and the scaling parameter, ao, respectively. The 

temperature range of the scaled data for each film is from T=32.4 K to 35.0 K (a;=0.15), from T=36.0 K to 37.6 K (a;=0.16), 
from T=28.2 K to 31.0 K (x-=0.18), from T=27.0 K to 29.0 K (a-=0.19), and from T=33.2 to 36.5 K (a;=0.20), respectively. 
Solid (dashed) lines are the 2D (3D) Gaussian scaling functions. 
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FIG. 8: Scaled data of (j)^ for the 2;=0.17 films with dif- 
ferent thicknesses. They are plotted in a temperature range 
from T=29.0 K to 30.8 K for the i=115 nm film, and from 
r=35.0 K to 36.6 K for the t=140 nm film. Solid (dashed) 
lines are the 2D (3D) Gaussian scaling function. Insets: Tem- 
perature dependence of dc resistivity of the corresponding 
films. 



FIG. 6: Temperature dependence of the obtained scaling 
parameters, uoo-q (top), ujo (middle), and ctq (bottom). Left 
(Right) panels are for the optimally doped (the overdoped) 
LSCO. Solid lines in the left panels are calculations given by 
Eqs. (2) and (3) with u=0.67, z=2, d=3, suggesting the 3D- 
XY critical charge dynamics. Dashed lines in the right panels 
are similar calculations with d=2 and uz=1.5. 
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FIG. 7: Plots of the highest (or the lowest) temperatures 
where a critical slowing down in ujo (T) was observed as a func- 
tion of the Sr concentration (hatched area). Difference in sym- 
bols represents the difi'erence in the observed critical behavior, 
which is classified into three groups of 2D-Xy (squares), 3D- 
XY (circles), and 2D- "U" (triangles) . Below the dashed line, 
the Nernst signal was observed.— 




FIG. 9: Scaled data of (j>a for the a;=0.18 films with (a) 60 nm 
thick and (b) 240 nm thick. They are plotted in a temperature 
range from r=28.2 K to 31.0 K for the f=60 nm films, and 
from r= 35.0 K to 37.0 K for the f=240 nm films. Insets: 
Temperature dependence of dc resistivity of the corresponding 
films. 




FIG. 10: Dashed (solid) thick line represents a critical di- 
vergence of ^(T) = io{T/Tc - 1)°-®^ in the ZD-XY critical 

charge dynamics with =3.5 nm (20 nm). Dashed and solid 
arrows show temperature ranges where the critical behavior 
was observed for the two LSCO (a;=0.16 and a;=0.18, respec- 
tively). Thin solid lines show the thicknesses, t, of the two 
films. Note that the SD-XY critical charge dynamics of the 
a;=0.16 film was not cut off by the film thickness (t=140 nm) 
in the observed temperature range. This suggests that is 
smaller than 3.5 nm. On the other hand, if the 2D critical be- 
havior of the a;=0.18 film in the observed temperature range 
was attributed to the ?>T>-XY critical behavior, which was cut 
off by the film thickness (t=240 nm), must be larger than 
20 nm. 
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FIG. 12: Temperature dependence of (a) ujo, (b) ctq, and (c) £, oc l/y/Uo for the x = 0.07 film. Solid lines in each panel are 
calculations given by Eq. ((9} with ^o=3 nm. 
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FIG. 13: Scaled data of (j)a (upper panels) and \a\ (lower panes) for the a;=0.16 film in zero and finite magnetic fields. 




FIG. 14: Temperature dependences of the obtained scaling parameters, (a) Wofo, (b) wo, and (c) ctq, for the x=Q.lQ film in 
zero and finite magnetic fields. Solid lines are calculations given by Eqs. (2) and (3) with f=0.67, 2=2, and d=3, suggesting 
the ZD-XY critical charge dynamics. 
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FIG. 15: Scaled data of (j)a (upper panels) and |cr| (lower panes) for the a;=0.18 film in zero and finite magnetic field. 




FIG. 16: Temperature dependences of the obtained scaling parameters, (a) uioao, (b) wo, and (c) cro, for the a;=0.18 film in 
zero and finite magnetic fields. Solid lines are calculations given by Eqs. (2) and (3) with i'z=l.h and d=2. 
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FIG. 17: Plots ol the applied magnetic field for the optimally 
doped LSCO (a;=0.16) and the overdoped LSCO (x=0.18) as 
a function of 1 — Tc{B)/Tc{Q). Sofid and dashed lines are 
calculations given by Eq. (llOp with v=Q.&7 and ^'=0.9, re- 
spectively. Inset: Plots of the applied magnetic field, B, as 
a function of the critical temperature, Tc{B) for each film, 
which was determined through the dynamic scaling analysis 
at each magnetic field. If vortex dynamics can be neglected, 
this plot agrees with the temperature dependence of Bc2(T). 
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FIG. 18: Scaled data ol (pa (upper panels) and \a\ (lower panels) for the x=Q.n film in zero and finite magnetic fields. Solid 
(dashed) lines are the 2D (3D) Gaussian scaUng functions. 




FIG. 19: Temperature dependences of the obtained scaling parameters, (a) uioao, (b) wo, and (c) cro, for the x=Q.17 film in 
zero and finite magnetic fields. Solid (dashed) lines are calculations given by Eqs. (2) and (3) with vz=l.h and d=2 (i/=0.67, 
2:=2, and d=3). 
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(a) independent vortex pair (b) vortex ring 

FIG. 20: Schematic drawings of (a) independent vortex- 
antivortex pair (b) vortex ring in a layered superconductor. 




FIG. 21: Schematic phase diagram of the 2D quantum an- 
tiferromaenet. Lines fal. Cb) and (c) renresent fo in each 
regioi 
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Carrier density 



FIG. 22: Schematic phase diagram of LSCO together with 
our resuhs about the critical charge dynamics. Solid circles 
are the phenomenologically well-established QCPs. An open 
circle near the optimal doping is a hidden QCP assumed to 
explain our findings. The region between two dashed lines is 
the quantum critical region. Thick solid line at T=Q is a long 
range ordered region. 
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FIG. 23: An example of the breakdown of the dynamic scal- 
ing for the (j>cr of the thicker film (t=120 nm) with x—0.18. 
Inset; Temperature dependence of do resistivity of the film 
compared with that of the thinner film {t=60 nm). 




FIG. 24: Temperature dependence of 1/r, where r is the 
correlation time. Solid symbols are the data of LSCO films 
{x=0.07, 0.16, 0.18) in this work, and open symbols are the 
data of NbN films {t=50 nm, 150 nm) in Ref. A solid 
straight line is I/tgl given by Eq. (I12|l . 



